Cross sections for short pulse single and double ionization of helium 
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In a previous publication, procedures were proposed for unambiguously extracting amplitudes for 
single and double ionization from a time-dependent wavepacket by effectively propagating for an 
infinite time following a radiation pulse. Here we demonstrate the accuracy and utility of those 
methods for describing two-photon single and one-photon double ionization of helium. In particular 
it is shown how narrow features corresponding to autoionizing states are easily resolved with these 
methods. 
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I. INTRODUCTION 

The development of attosecond-pulse radiation sources 
[H offers the prospect of a new class of pump/probe ex- 
periments that can in principle explore the effects of elec- 
tron correlation in atoms and molecules on ultrashort 
time scales [l, 3 3- The analysis and interpretation 
of such experiments will necessarily involve state-of-the- 
art time-dependent, non-perturbative theoretical meth- 
ods and advanced supercomputing resources. For exam- 
ple, if the pump step excites autoionizing states of the 
target whose lifetimes are much longer than the pulse du- 
ration, then traditional time-dependent methods might 
require prohibitively long integration times to compute 
meaningful ionization probabilities. 

In a previous paper Q, hereafter referred to as paper 
I, we outlined a procedure for extracting the amplitudes 
for ejecting electrons of particular energies and directions 
from a quantum wavepacket at the end of a short pulse, 
while the electrons are still interacting with the target 
nucleus and each other. The basic idea was to solve the 
time-dependent Schrodinger equation over the finite pe- 
riod of time when the pulse was acting on the target and 
to then effectively propagate the solution to infinite time 
by using the propagated wavepacket as the source term 
in a time-independent driven Schrodinger equation with 
the field-free Hamiltonian. The method proposed in I 
was illustrated with computation of one-and two-photon 
ionization cross sections for atomic hydrogen and we out- 
lined the theory for extending the method to two-electron 
targets. 

In this paper, we demonstrate the viability of the 
approach with computations on atomic helium. We 
will show that the present method allows us to ex- 
tract fully differential ionization probabilities over the 
entire bandwidth of the pulse and to resolve structures 
arising from relatively long-lived autoionizing states, 
which might require prohibitively long propagation times 
with traditional time-dependent approaches. Although 
the methodology can be applied with arbitrary field 
strengths, we confine our attention here to low inten- 
sity fields so that we can compare the present method 



with the results of other studies that calculated one- and 
two-photon ionization cross sections in the perturbative 
limit. 

The outline of this paper is as follows. The theory is 
outlined in Sec. [TTl beginning with a derivation of the 
driven equation and the extraction of ionization ampli- 
tudes, followed by explicit formulas for one- and two- 
photon cross sections for single and double ionization and 
a brief description of the computational procedures we 
employ. In Sec. IIIII wc present results for one- and two- 
photon single ionization of helium, while Scc. lIVI prcsents 
results for one-photon double ionization. We conclude 
with a brief discussion. 



II. THEORY AND COMPUTATION 

A. Time-dependent Schrodinger equation and 
extraction of ionization amplitudes 

The methodology we use is fully detailed in I and 
so only the essentials are repeated here. We assume the 
atom, initially in its ground state, is subjected to a time- 
varying pulse that starts at i = and ends at t = tfinai- 
To track the time evolution of the wave function dur- 
ing this period, we solve the time-dependent Schrodinger 
equation: 



(1) 



At t = tfinai, the time- varying field ends and the wave 
function continues to evolve under the time-independent 
atomic Hamiltonian; this time evolution can be written 
explicitly as: 



(2) 
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We next define a scattered wave ^'^c by taking the Fourier 
transform, from tfinai to infinity, of Eq. Q 
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or, equivalently, 
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Thus the scattered wave, from which we will extract 
physical information, satisfies a driven Schrodinger equa- 
tion in which the propagated wavepacket at the end of 
the pulse appears as the source term. 

To construct amplitudes for single and double ioniza- 
tion, we begin by formally expanding the wavepacket, at 
the end of the pulse, in the complete set of bound, single 
and double continuum eigenstates of the target Hamilto- 
nian: 

5'(ri,r2,<final) = 

-'/'boundCri, r2) -I- V'singlo(ri, r2) + V'd oublc 
= '0bound(ri,r2) 



dfc3C(k„)Vk (ri:r2) 



(5) 



+ j dklj &fC(ki,k2)Vk„k.(ri,r2), 

where ■!/'bound(ri, i"2) contains the contributions from the 
bound states of the target, n runs over the bound states 
of He+ and the coefficients C(k„) and C(ki,k2) are am- 
plitudes for single and double ionization, respectively. As 
we showed in I, substituting Eq. ([5]) into the expression 
^'sc = G+\['(tfinai), and using the asymptotic form of the 
full Green's function, allows us to write the asymptotic 
forms of ^sc for single [1] and double ionization [J: 



J{kri+{Z/k} ln2fcri) 
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where 
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(j>n is a bound state of Hc+ and p ~ + ^2 • 

By solving the driven equation ^ using exterior com- 
plex scaling (ECS), pure outgoing boundary conditions 



are automatically imposed on the scattered wave func- 
tion. And having identified the amplitudes C(k„) and 
C(ki,k2) in the asymptotic form of the scattered wave, 
their explicit evaluation is done in terms of the following 
surface integrals @: 



C(k„) = i J [0^*(ri)C(r2)(V*sc(ri,r2)) 

-vI/,,(ri,r2)V(0k*(ri)C(r2)] -d^, 
for single ionization, and 



(8) 



C(ki,k2) 



-vI/,c(ri,r2)V(0k;(ri)0-;(r2))] -rfS, 
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for double ionization, where V = (Vi,V2) and x is a 
volume-dependent phase that makes no contribution to 
any physical observable The testing functions cp^ 
arc momentum-normalized Coulomb functions with a nu- 
clear charge Z = 2. 



B. One-photon cross sections 

The amplitudes for ionization extracted via Eq. ([8|) 
or Eq. ^ will generally depend on the parameters (in- 
tensity, bandwidth, etc.) of the radiation pulse that 
produced the wavepacket being analyzed. However, if 
the intensities are such that time-dependent perturba- 
tion theory gives an accurate description of the physical 
process, then the amplitudes can be used to construct 
one-photon cross sections and, if the pulse durations are 
not too short, two-photon cross sections, over the range 
of energies within the bandwidth of the pulse. 

In the dipole approximation, the laser-atom interaction 
in the velocity gauge is given in terms of the electron's 
momentum operator p and the vector potential A by 
U{r,t) = (e/mc)A(i) • p. For a photon energy ui and a 
total pulse duration T, A(t) may be written 



A{t) 



AoFUt)e te [0,T] 
elsewhere 



(10) 



where e is the polarization vector. We choose a sine 
squared envelope for the time and frequency dependence 
of the pulse, F^{t), 



sm~ I —t J sin((jjt) , 



(11) 



One-photon ionization cross sections are obtained 
using first-order time-dependent perturbation theory 
(TDPT). In first-order TDPT, the transition amplitude, 
C^^ , between an initial state of energy Ei and a con- 
tinuum final state of energy Ef, caused by a pulse of 
duration T characterized as in Eqs. pOl) and pT|) is 



e-p|$,)Fi"(t^,Ai;,T) 



(12) 
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where AE ={Ef ~ Ei) and 



C. Two-photon cross sections 



Jo 



(13) 



(r2(w - AE/hY ~ 47r2)(tj - AE/h) 



and ($J |e • p|(f>i) is the dipole matrix element connecting 
the initial and final states. 

The cross section for one-photon single ionization, leav- 
ing the ion in state n is 



I(*kje-Pl<i'« 



(14) 



dn mAE 
while the one-photon double ionization cross section is 

Combining Eqs. ([14]) and ([15]) with Eq. ([12]) then gives 
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and 
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aAE \Ao\^\F^^{uj,AE,T)\^' 
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4Tr^kik2e'^m'^ 



|C(ki,k2)|^ 



dEidnidn2 ah^AE \Ao\^\F^'^{uj,AE,T)\'^' 

. . 

Thus the factorability of the transition probability in 
first-order TDPT allows us to extract the one-photon dif- 
ferential single and double ionization cross sections from 
a single pulse within its bandwidth, defined as the range 
of photon frequencies for which F^^ {to, AE,T) is appre- 
ciably nonzero. 

It is worth pointing out that for long pulse durations 
\F^'^{uj,AE,T)\'^ tends to {3nT/16)6{uj ~ AE/h) 
This identity may be used to define cross sections at 
uj = AE/h via Fermi's golden rule, but this is only mean- 
ingful when the long time limit has been reached. Equa- 
tions (fTB|) and (fT7|) . on the other hand, properly define 
the cross sections for any length pulse - a point that has 
evidently been missed in the recent literature p^ . 



For a two-photon absorption process, we can use 
TDPT in second order to write the amplitude for a tran- 
sition between an initial state of energy Ei and a contin- 
uum final state of energy Ef as 
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X F^^{Ef,E^,E,,u;,T), 



(18) 



where the sum m is over all the eigenstates of the target. 
The coefficient F^'^{Ef,E^,E„uj,T) is given by 



F^^iEf,E^,E,,u,T) = 

- f dt'e^^'^^/'" sin2(t'7r/T) 



X 1 /" (it"e*(^^^""/'*— )*"sin2(t'V/r) 
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where AEij — Ei — Ej 



To connect Eq. (jlSp with the familiar expression for the 
two-photon cross section, we define a reduced coefficient 
or " shape function" : 



F^'^{Ef,Ern,E„LO,T) 



{E, + AEf,/2 - Ern)F^'^{Ef, E^, E,,iu, T). 



(19) 



As we explained in I, we have found that if the photon 
frequency lj is not too close to being in resonance with a 
transition to one of the intermediate states, m, then the 
shape function F'^'^ is well approximated by an expression 
that does not depend on the energies of the intermediate 
states in the sum in Eq. (jlSp and which becomes exact 
in the long T limit: 



F'^^{Ef, E,n, E,,uj, T) « :S{Ef, E,, T) 



(20) 



(2w - AEf,) [T^{2uj - AEf,Y - 20tt^T^{2uj - AEfif + Utt^] 



giving We can now connect the amplitudes we extract from 

the wavcpacket using Eqs. (I8|) or (l9|) with differential 
iuAqS. 1 cross sections for two-photon ionization. The two-photon 



[E, + AS/,/2 - H) 
X ^{Ej,E,,uj,T) 



(21) 
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single ionization cross section is given by the expression, 



(22) 



while the two-photon double ionization cross section is 
given by 



da 



2uj 



(27r)3fcifc2to2 
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(23) 



Combining the previous two cross section definitions with 
Eq. (ini) then gives 



da^ 



|C(k„) 



dn 
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and 



(24) 



da^ 



(2Trfhkik2a^ 



|C(ki,k2) 



dE,dnM iAE,j2ym^ mEf,E.,.,TW 

(25) 

Once again, the factorability of the transition probabil- 
ity allows us to extract cross sections from a single pulse 
within its bandwidth, but we must emphasize that, in 
contrast to the one-photon case, this factorability is ap- 
proximate and, as we shall see, breaks down when the 
pulse width is very short or may require very long prop- 
agation times to resolve the energy dependence of the 
cross section when the photon frequency is close to being 
in resonance with a discrete intermediate state. 



D. Implementation 

We employ the same computational techniques here 
that we have used our recent work on two-electron prob- 
lems [ll|, [13 • The two-electron wave function is first 
expanded in products of spherical harmonics 



*(i-i,i-2,t) 



EE 



1 

^I'^s (26) 

imi m2 v-2 ) • 



Wc include all /m-pair configurations that can be con- 
structed using some specified value of Imax- Integra- 
tion over the angular variables then gives a set of cou- 
pled equations for the two-dimensional radial functions 
iphmi ,i2m2 {fi I ''2 , ■ The radial degrees of freedom are in 
turn discretized using a finite-element, discrete variable 
representation (FEM-DVR) with a product basis of Lo- 
batto shape-functions [l^. The value of Imax, sls well as 



the size and density of grid points required for conver- 
gence, will depend on the photon energy as well as the 
process under consideration. 

Exterior complex scaling of the radial coordinates. 



Ro + {r- Ro)e 



i9 



r<Ro 
r > Ro, 



(27) 



defines a radius Rq beyond which the radial coordi- 
nates are complex-valued. The round-state wave func- 
tion ^pQ, which serves as the initial wavepacket for the 
time-propagation, is obtained by diagonalizing the field- 
free Hamiltonian on a relatively small portion of the real 
grid (rmax ~ 30 Bohr), with configurations representing 
total angular momentum L = 0. This wavepacket is then 
propagated on the real part of the grid over the duration 
of the pulse. The time-propagation is carried out using 
a Cranck-Nicholson propagator, 

(1 - ^i7)*(ri,r2,i + At) = (1 + '-^H)^{r,,r2,t) 

(28) 

with a time step At on the order of 10 ^ atomic units. 
Since the Hamiltonian is time-dependent, Eq. ((28)) re- 
quires inverting a large matrix at each time step. How- 
ever, since the wavepacket changes little with each At, we 
have found that an iterative solution at each step, which 
requires only matrix-vector multiplications, converges in 
several iterations. Since the time-propagation is carried 
out on the real part of the grid, Rq must be chosen large 
enough to contain the spreading wavepacket over the du- 
ration of the pulse and avoid unphysical reflections from 
the grid boundaries. The value of Rq required will gen- 
erally increase with photon energy; in these calculations, 
we found that Rq ~ 130 Bohr was sufRcient for the range 
of photon energies considered. 

The wavepacket at the end of the pulse serves as the 
driving term for the scattered wave equation (Eq. Q). 
This equation is solved on the full exterior scaled grid 
and provides the scattered wave from which the physical 
amplitudes are extracted, as outlined above. We reiter- 
ate that the time propagation is carried out only once 
for a particular laser pulse and then the scattered wave 
equation can be solved for any energy within the pulse 
bandwidth. 



III. SINGLE IONIZATION OF HELIUM BY 
ONE OR TWO PHOTONS 

We restrict our calculations to weak fields where per- 
turbation theory can be applied and, therefore, a cross 
section defined. In this way, we can check the accuracy 
of the method by comparing with existing calculations 
and experimental measurements. 

All results reported here were obtained with an inten- 
sity 1= 10^^ W cm~^, which is high enough to provide 
relatively large ionization rates for one and two photon 
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FIG. 1: (Color online) Energetics of helium ionization. 
Schematic representation of (a) one-photon single ionization 
process, (b) and (c) two-photon single ionization processes; 
and (d) one-photon double ionization process. 



transitions and low enough to keep the processes within 
the perturbative regime ■ 

Firstly, we look at single ionization by absorption of 
one or two photons using pulses of duration T=0.9 fs. 
Over the range of photon energies considered, calcula- 
tions are converged by including pairs of spherical har- 
monics built with individual angular momenta up to 
Imax — 3 and leading to total angular momenta up to 
L = 2. In order to be consistent, we construct our ground 
state with the same maximum value for individual angu- 
lar momenta. In the weak field limit, optical selection 
rules apply, so only channels with total symmetry ^S'', 
and (i = 0, 1 and 2) are accessible. 



A. Resolving autoionization resonances in 
one-photon single ionization 

Figure [T] shows a schematic representation of one- and 
two-photon ionization processes in atomic helium. One- 
photon ionization, labeled (a), is possible with photon 
energies above 24.6 eV. For two-photon single ionization, 
with a threshold at half this energy, we can distinguish 
above-threshold ionization (ATI), labeled (b), in which 
the first absorbed photon is above 24.6 eV, from case 
(c), where two photons are required to ionize the target. 
For photon energies above 79.01 eV, one photon double 
ionization, labeled (d), is possible. 

The total cross section for one-photon single ionization 
is plotted in the top panel Fig. [2] as a function of photon 




Photon energy (eV) 

FIG. 2: (Color online) Top panel: One-photon single ioniza- 
tion cross sections, in units of megabarns, as a function of 
photon energy. lMb=10~^*cm^. Thin black full line: time 
independent perturbation theory results from reference [isl ]. 
Thick lines: Present results for T= 0.9 fs, / = 10^^ W cm^ and 
different number of cycles. Inset in top panel: Enlargement in 
the region where doubly excited states are resolved. Bottom 
panel: Squared amplitudes, in atomic units, from which the 
total cross sections in the top panel were extracted. 



energy. These cross sections were extracted from calcu- 
lations using different central frequencies of the field, i.e. 
different numbers of optical cycles. For a pulse duration 
T= 0.9 fs, 6 cycles correspond to a pulse of central fre- 
quency around 27.6 eV, 8 cycles correspond to 36.76 eV, 
and so on. For the cross sections depicted in Fig. ^ the 
corresponding squared amplitudes, appearing in the bot- 
tom panel, capture the energy bandwidth of the pulse, 
i.e., the Fourier transform of F^{t) defined in Ea. pTI) . 

The results plotted in Fig. [5] are practically indis- 
tinguishable from time independent perturbation theory 
calculations (black thin line) and experimental measure- 
ments given in ref. (isj . Agreement is excellent even in 
the region between 60 and 65 eV where one-photon ab- 
sorption can populate doubly-excited states of sym- 
metry, the first four of which have been labeled in the 
inset of the top panel. Cross sections in this region were 
obtained from calculations using pulses of 0.9 fs with 12 
and 14 optical cycles. 

We should point out that conventional time-dependent 
treatments would require long propagation times to fully 



6 



resolve the autoionizing structures seen in Fig. [2]- greater 
than 6ps, for example, in the case of the resonance 
near 62.7 eV. With the present method, we can obtain 
such results with much shorter propagation times be- 
cause the time span from the end of the pulse to infinity 
- when the two electrons are still interacting - is handled 
exactly by Eq. SI 



B. Two photon single ionization 

We have also explored single ionization of helium by 
two-photon absorption. Total cross sections for this pro- 
cess are plotted in Fig. [3] as a function of photon energy. 
Given the optical selection rules, only states with and 
"'^D*' final symmetries will be populated from the ground 
state. The corresponding amplitudes appear in the bot- 
ton panels for various pulses with different optical cycles. 

Structures appearing in the total cross section for pho- 
ton energies above 28 eV correspond to doubly excited 
states of symmetries ^S'^ and "'^D'' (labeled in the figure) 
which decay to the continuum after the pulse is turned 
off, as we have discussed in the previous subsection. 
We note that, at these photon energies, above thresh- 
old ionization (ATI) processes are taking place [sec en- 
ergy scheme of Fig. [1] process labeled (b)]. The positions 
and widths of these autoionizing states are in reasonable 
agreement results obtained from time- independent per- 
turbation theory calculations in refs. [l6l. ITtI. Iist . 

We must point out that our calculation with a 0.9 fs 
pulse does not reproduce any of the structures between 20 
and 24.5 eV found with time-independent perturbation 
theory [l^. At these energies, bound excited states of 
helium are populated by one-photon absorption [process 
represented in Fig. [T] as (c)], which leads to divergences 
in the total cross section when the pulse length goes to 
infinity. Since we are using a finite length pulse with 
an energy bandwidth wider than the energy spacing be- 
tween these intermediate bound states, these cannot be 
resolved. 

Figure 0] shows two-photon single ionization cross sec- 
tions extracted from calculations using pulses of different 
durations. These results highlight the completely differ- 
ent nature of the structures arising from single-photon 
absorption by intermediate bound states and the struc- 
ture associated with doubly excited states. The latter 
appear at photon energies above 25 eV where the cross 
sections are seen to be invariant to increasing the pulse 
length. The doubly excited states are first populated, in 
the presence of the field, by two-photon absorption and 
decay later in time. By contrast, the track of intermedi- 
ate states between 19 eV and the first ionization thresh- 
old only begins to appear in the extracted cross sections 
for long pulse durations. With a short pulse duration of 
0.45 fs, the effective bandwidth is Aw = 27r/T= 9.2 eV, 
which is too broad to resolve any structure below 25 eV. 
For a 3 fs pulse, the bandwidth is ~ 1.4 eV and we see 
structure beginning to develop below 25 eV. Much longer 



pulse durations would be required to fully resolve these 
structures. 
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FIG. 3: (Color online) Upper panel: S and D symmetry com- 
ponents of the total cross section for two-photon single ion- 
ization of helium, calculated with pulses of duration of T=0.9 
fs and / = lO^'^ W cm~^ and varying number of optical cycles 
(4, 5, 6 and 7). Lower panel: Squared amplitudes, in atomic 
units, from which the total cross sections in the top panel 
were extracted. 



So to conclude this section, we reiterate that in the 
two-photon case, there are two different issues with re- 
spect the pulse duration that must be carefully consid- 
ered. First, as we have just discussed, is the fact that 
with finite pulses intermediate bound states will be re- 
solved only if the energy bandwidth of the pulse is nar- 
rower than the structures in question. This is simply the 
physics of the problem. The second point is that the 
factorability of a "shape function" from the transition 
probability for two-photon absorption relies on an ap- 
proximation (Eq. (j20|l ) that becomes unreliable for very 
short pulses. Indeed, careful examination of Fig. 3] shows 
that in the case of a 0.45 fs pulse, there are small errors in 
the computed cross sections even above 25 eV. We hasten 
to remark that these issues are only relevant when one 
wishes to compute cross sections to compare with the re- 
sults of time-independent perturbation theory. Whether 
or not we operate in the perturbative limit, the ampli- 
tudes computed for any pulse length or field intensity will 
still give correct ionization probabilities. 
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FIG. 4: (Color online) Dependence of extracted two-photon 
single ionization cross section on pulse length. 



IV. ONE PHOTON DOUBLE IONIZATION 

A key test of the method is the extraction of total and 
differential cross sections for double ionization processes. 
In Fig. [3] we show the total cross section for one-photon 
double ionization (process labeled (d) in Fig. [5]) . The to- 
tal cross section was calculated for a wide range of photon 
energies by using only two different wavepackets propa- 
gated with a 250 attosecond pulse. The agreement with 
experimental results by Samson [l9| is excellent. Once 
again, it is the exact factorability of the amplitude in the 
one-photon case that allows us obtain accurate results 
with such a short pulse duration. 

Since one-photon double ionization processes are ex- 
tremely sensitive to electron correlation effects in both 
initial- and final-states , single and triple differential 
cross sections can be sensitive to higher values of elec- 
tron orbital angular momentum. We have checked that, 
for the range of energies considered here, convergence is 
achieved by including Imax = 4 for individual angular 
momenta in the spherical harmonic basis for the propa- 
gating wavcpackct as well as for the ground state. 

In Fig. [HI we plot the single differential cross section 
(SDCS), as a function of energy sharing, for four different 
total energies. Our results are in excellent agreement 
with previous experimental measurements |20ll. as well 
as the calculations of Colgan and Pindzola [21[ . 

In Fig. [71 we show triple differential cross sections 
(TDCS) at 20 eV above the double ionization threshold, 
where the total cross section reaches its maximum value. 
Most of the experiments, including those shown here, 
have been performed in coplanar geometry, i.e., with the 
polarization vector and both momenta ki and k2 in the 
same plane. For the TDCS, the time propagation is car- 
ried out with a 6 cycles pulse (central frequency of 99.98 




Excess Photon Energy (eV) 

FIG. 5: (Color online) Upper panel: total cross section, in 
units of kilobarns, for one-photon double ionization of he- 
lium. 1 kb=10~^'^cm'^. Solid curve: present results; dots: 
experimental results of Samson [l^. Lower panel: squared 
amplitudes, in atomic units, from which the total cross sec- 
tions in the top panel were extracted. 




E,/E 

FIG. 6: (Color online) Single differential cross sections, in 
units of kilobarns per eV, for one-photon double ionization of 
helium at different total energies for the ejected electrons. 



cV), and the extraction for a total absorbed energy of 
99.01 cV. Results are shown for three different fixed val- 
ues of ejection angle for electron 1(30°, 60° and 90°) and 
three different energy sharings. We find excellent agree- 
ment with the absolute experimental results of Brauning 
et al. [2^, as well as the convergent close-coupling (CCC) 
calculations of Kheifets and Bray [l^l . 

Figure [8] shows TDCS results at 60 eV above thresh- 
old for equal energy sharing of the ejected electrons. The 
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FIG. 7; (Color online) Triple differential cross sections 
(TDCS), in units of barns per eV per steradian^, for one pho- 
ton double ionization of helium at 20 eV above threshold. 1 
barn=10~2''m^ Solid curves: present results; broken curves; 
CCC results of Kheifets and Bray [S^l; points: experimental 
results of Brauning et al. [2^ . 



experimental results in this case [2J] were internormal- 
ized but not absolute, so we normalized them to our 
present results for comparison. The CCC theoretical 
results for this case arc seen to agree with our calcula- 
tions in shape, but are ~20 percent smaller in mag nitude. 
Time-dependent close coupling calculations [2l| are in- 
distinguishable from our results (for this reason these re- 
sults are not plotted). 

Finally, we present TDCS results for a photon energy 
of 85 eV, which is only 6 eV above the ionization thresh- 
old. The TDCS are plotted in Fig. [H aloiig with ab- 
solute experimental results of Dorner et al. [25[ . In the 
experiment, the data was binned over finite ranges of en- 
ergy sharing( 0-0.2, 0.5 and 0.8-1) and polar angle 6i 
(40-65°) for the fixed electron. Our calculations are for 
fixed energy sharings of 0.1, 0.5 and 0.9 and we fixed the 
polar angle 9i at 45°. In Fig. [51 each row corresponds 
to a different energy sharing while the columns corre- 
spond to a different value of azimuthal angle (j) between 
the fixed and varied electrons. The calculations were 
done for (/) = 0°, 30° and 60°, while in the experiment (j) 
was binned (0-20°, 40-65° and 45-90°). We find excellent 
agreement with experiment. Also shown in Fig. [5] are the 
results from ref. (25[ obtained using fourth-order Wannier 
theory [2^ , which do not provide absolute cross sections 
but evidently do reproduce the shapes of the TDCS at 
these low energies. 
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FIG. 8: (Color online) As in Fig. [3 at 60 eV above threshold. 
Solid curves: present results. Dots: experimental results [24| . 
normalized to the present results. Broken curves: CCC from 
ref. [Hi 
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FIG. 9: (Color onhne) As in Fig. [3 at 6 eV above threshold. 
Solid curves: present results; experimental results (dots) and 
fourth-order Wannier theory results (broken curves) are from 
ref. [H. 



V. DISCUSSION 

We have demonstrated that the method we proposed in 
I can indeed be applied to the study of single and double 
ionization processes in two-electron systems. By working 
with relatively low field intensities where perturbation 
theory is expected to be valid, we have been able to show 
that the cross sections we calculate for atomic helium 
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are consistent with earlier results obtained using time- 
independent methods, as well as experimental data. By 
propagating a wavepacket in the presence of a single pulse 
over its time duration, we can extract cross sections over 
the entire bandwidth of the pulse, even in the ATI region. 

Having demonstrated that we can reproduce the re- 
sults of time-independent calculations with this method, 
we hasten to add that our main purpose in develop- 
ing this approach is to provide an efficient and reliable 
method for the exploration of problems that are difhcult 
or impossible to study with time- independent techniques. 
Such problems include studies of above-threshold, two- 
photon double ionization, which are far from straightfor- 
ward with time- independent perturbative methods ^T^ . 
short-pulse, intense field studies where perturbation the- 



ory is not valid, and simulations of two-color, pump- 
probe experiments which require an approach that is ex- 
plicitly time-dependent. 
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